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Abstract 

We propose a method for finding gaps in the spectrum of a differen- 
tial operator. When applied to the one-dimensional Hamiltonian of the 
quartic oscillator, a simple algebraic algorithm is proposed that, step 
by step, separates with a remarkable precision all the energies even for 
a double-well configuration in a tunnelling regime. Our strategy may 
be refined and generalised to a large class of ld-problems. 

PACS numbers: 03.65.Db, 02.30.Tb, 02.60.Gf, 02.60.Pn 

1 Introduction 

Obtaining some physically interesting piece of information on a spectrum 
of a differential operator is generically a difficult problem. For systems in 
low dimension, say one or two, though numerical approaches can usually 
compute the spectrum to a high precision, very few is usually known when 
demanding the standard criterions of rigour of mathematical physics. 

The aim of the present article is to propose a method for finding gaps 
in the spectrum of a differential operator. We will mainly focus on the sta- 
tionary one-dimensional non-magnetic Schrodinger equation with a smooth 
potential V{x) defined on the whole real axis, 

-ig + V„-*, <1) 

but it can be easily understood that our ideas, which will be explained in the 
first part of this article (§ [2]), can be extended to a wider family of ordinary 
differential equations. For the sake of definiteness, in the second part (§ [3]) 
we will illustrate more concretely our method on the quartic oscillator (jlip 
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where the numerical "exact" spectrum will serve us as a touchstone. We will 
show that we are able to obtain gaps in the energy spectrum, i.e. intervals 
where there is no eigenenergy, with purely algebraic manipulations involving 
no integrals. An algorithm is explicitly constructed in the case of the quartic 
oscillator — but it remains robust under any kind of smooth perturbations 
- and provides more and more non-trivial gaps. As far as we know, in the 
huge mathematical physics literature (Bender and Orszag, 1978; Voros, 1983, 
for instance and their references) on the quartic oscillator no such gaps have 
been found. We will present strong evidence that, for still unknown reasons, 
this algorithm allows to resolve all the eigenenergies with a surprising pre- 
cision even though, remaining stuck to the simplest means, we are far from 
having exhausted all the possibilities that our method offers. Actually, some 
very subtle information on the spectrum can be gained since we can find a 
gap in between the tunnelling doublets, i.e. when the energies are below the 
height of the intermediate barrier in a double-well configuration. 

2 Guideline of the method 

2.1 General strategy 

The key idea of our method is to construct, for a given integer N ^ 1, a 
real function Jjsr(tp' ', <p, x, E) from a possible real square-integrable solution <p 
of ([1]) with energy E (the prime stands for the derivative d/dx) such that 



The real function Fn of both the spatial coordinate x and the energy E is ip- 
independent and is constructed from the potential V and its derivatives. An 
example of Jn is given by equation (0) below and some F^'s are given by §5\ 
-[8]). Condition (i) is the cornerstone of our method and, before justifying 
how it can be obtained (we will see that condition (ii) is not so restrictive), 
let us first explain how gaps in the spectrum may be obtained. 

When the conditions (i) and (ii) are simultaneously fulfilled, an im- 
mediate consequence is that the integral f^°° (Wx)) Fn(x,E)<1x van- 
ishes. This implies that, if E is truly an eigenenergy, the function x i— > 
Up(x)} Fn(x,E) should change its sign. If N is even, we obtain a ip- 
independent condition: for any fixed energy x \— > Fn(x,E) must change its 
sign on the real axis. We will see that for such a one-dimensional prob- 
lem, and for a given N, we still can choose Fn in a wide continuous set of 
smooth functions on the real axis. A forbidden value of E (i.e. E cannot be 
an eigenenergy) is obtained if we are able to construct a Fn that remains 
positive on the whole x-axis. Once this property is achieved, it remains 




(j N (tp'(x),(p(x),x,E)^ = (tp(x)) N F N (x,E) 
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stable under small perturbations within the set of -Fat's, and we obtain a 
whole interval where no eigenenergy can exist. 

When N is odd, a priori some non-trivial information can be extracted 
for the ground state only, if there is any. Indeed, it is known that its wave- 
function can be chosen positive for any Schrodinger equation of the form ([T]) 
(see for instance (Reed and Simon, 1978, XIII. 12)) whereas all the excited 
eigenfunctions do change their sign. No gap can be found but one can use 
this strategy to find upper and lower bounds on the ground state energy 
even in a multidimensional situation. The differential method presented in 
(Mouchet, 2005) corresponds to N = 1. 



2.2 Obtaining condition (i) and determination of F 1 



N 



In order to obtain condition (i), let us start with J being a homogeneous 
polynomial of degree N with respect to its first two variables (<//, </?): 



N 

. N—n i 



J N ((p'(x),ip(x),x,E) a n (x,E) (ip'(x)) n ((p{x)) n . (2) 



n=0 



The smooth functions {o- n } n £{o,...,N} wm be constructed in order to get 
condition (i) : Equation ([1]) allows us to eliminate ip" from the total derivative 
of Jat and then, the systematic cancellation of the coefficients of (tp*) n 
for N — n > leads to the relations 

1 TV — 77 + 1 

Vn€{0,...,JV-l}, a n+l = — a' n ^— 2 (V -E)a n ^, (3) 

n+l n + 1 

(to get a unified expression, we define a_i =0) with 

F N {x,E) =a' N {x,E) + 2(V(x)-E)a N ^{x,E) . (4) 

The recurrence relation (j3J) uniquely determines all the a n 's and Fjv from ao, 
which remains a free smooth function. For instance, we have 

F 1 = - 4 + 2(V -E)a , (5) 
F 2 =1 og' - 4(V — E)a' - 2V'a , (6) 

F 3 = -l4 v) + f(V-E)^ + fv'a' +(V"-6(V-Ef)a , (7) 

Fi =k ^ -\ {y - E) a " - \ v ' a ° + (-5 y " + f {v ~ E? ) a '° 

+ (-l v '" + Y v ' {v - E) ) a °- (8) 

The only condition on ao is that it must not increase too rapidly when \x\ — > 
+oo in order to get condition (ii). From the standard semiclassical analysis 



3 



(Messiah, 1991, for instance chap. VI), we know that if a bound state of V ex- 
ists, its wavefunction decreases exponentially as exp(— | f x y/2(V(x') — E) dx'\ 
when \x\ — ► +00; therefore the ansatz ([2]) will vanish at infinity as soon as 
I oto I becomes negligible compared to \ip\. The condition 



f x J2(v(x')-e) dx' 
\ao(x)\ <C e when \x\ — > +00 (9) 

is sufficient and not too demanding. 

We do not loose generality with the form ^§ . In fact any smooth function 
J [if 1 , (p, x, E) will inevitably lead to the hierarchy of functions {Fn} n>1 . To 
understand this, consider one monomial c(ip') n (f m , where c is a function of 
x, in the Taylor expansion of ^J((p' ,<f,x,E) with respect to (<p',(f) once 
tp" has been substituted by 2(V — E)(p. This monomial can be written as 

c (y/)> m = _ -J— (c^T- 1 )'^ 1 + — TTT~ Wr"V m+1 ) 
m + 1 m + ldx ' 

= — \d (ip') n ~ V m+1 + 2(n - 1)(V - E)c (^) n "V m+2 l 

m + 1 L J 

+^ T ^-(c(^r-v m+1 ). (10) 

m + 1 dx 

By using this type of identity an adequate number of times we can sys- 
tematically absorb all the powers of ip' in a total derivative while keeping 
the homogeneity in (<//,<£>). This procedure can be pursued until we ob- 
tain c(ip') n ip m of the form cip n+m + 4z(j((p',(p,x,E)) where c and j are 
smooth functions. A redefinition of J — » J — j for each monomial of de- 
gree n + m = N leads to condition (i) . Rather than the starting point ([2]) , 
we could have started from Jjv(y' ,<p,x, E) = ao((p') N and, working modulo 
a total derivative, repetitions of procedure (fTU|) would have led to condi- 
tion (i). For small N, it can been checked that the Fn thus obtained are 
also given by jS]-[S]). 



3 Application to the quartic oscillator and gener- 
alizations 

To illustrate the efficiency of our approach, in this section we will consider 
a quartic potential, which can always be reduced in appropriate units to 

V(x) ^ & ^ + \x\ (11) 

for s real. The associated energy spectrum is purely discrete and bounded 
from below by the minimum min(V). The energies can been computed 
numerically without any difficulty by diagonalizing the Hamiltonian in the 
standard basis of the eigenstates of an harmonic oscillator. They are given by 
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the continuous black lines in figured! For large positive s, the bottom of the 
spectrum {E n \ n = 0, 1, 2, . . . } tends to the harmonic spectrum (n+l/2)y / i. 
For s < 0, we get a double-well configuration whre the central barrier reaches 
its maximum at 0. When s decreases from 0, we observe the pairing of the 
energies (E2 n , -E^n+i) into doublets that characterize tunnelling from one 
well to the other. For s ^ sq — —2.0481 the two first states Eq and E\ are 
both negative, therefore below the energetic barrier, and 1/(E\ — Eq) rep- 
resents the tunnelling oscillation period (recall that we are working in units 
where h = 1) of a state initially localised in one well and constructed from a 
linear combination of the symmetric/antisymmetric eigenstates associated 
with Eq and E\ respectively. As — s increases, a standard semiclassical 
analysis shows that Ei — Eq behaves like 2 11 < /4 |s| 5 < /4 exp(— |2s| 3//2 /3)/y / 7r for 
large — s (Garg, 2000, § V, for instance). This exponentially small splitting 
is a highly non-trivial piece of spectral information to obtain by approximate 
methods. Nevertheless, we will show that for a given positive or negative s, 
we actually can exhibit a gap between Eq and E\. 

For simplicity, and in order to keep the computation tractable with ele- 
mentary algebraic manipulations, we will work with N = 2 and take ao of 
the form 

a (x) = P(x)e- Xx2 / 2 , (12) 

where A is real (not necessarily positive, see ([9])) and P is a real polynomial 
in x (to lighten the notations we leave the -E-dependence implicit. Note that 
P may depend on A as well). Condition Q is largely satisfied. From ex- 
pression Q, we have F2(x) = Q(x) exp(— Xx 2 /2) where Q is the polynomial 

Q( x ) = p"'( x ) - 3XxP"{x) - (2x A + (4s - 3A 2 )x 2 + 3A - 8E\p'(x) 

+ (2Ax 5 - (A 3 - 4sA + 4)x 3 + (3A 2 - 8EX - 4s)x) P(x), (13) 

that must change its sign if E is an eigenenergy. To find a forbidden value 
for the energy we construct an even Q and therefore start from an odd P: 

M 

P{x) = Y J P m x 2m+1 . (14) 

m=0 

Then 

M+3 

Q(x)=^2q m x 2m , (15) 

m=0 

where the coefficients {q m } are expressed in terms of the {p m } via the re- 
currence relation 

q m = (2m + 3) (2m + 2) (2m + l)p m+ i - (2m + l)(sX(2m + 1) - 8E^ Pm 

+ (2m(3X 2 - 4s) - SEX'jpm-i - (A 3 - 4sA + 4m - 2)p m _ 2 + 2Ap m _ 3 , 

(16) 
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for m € {0, . . . , M + 3} (we define p m = for m < or m > M). The 
simplest way to control the sign of Q is to reduce it to a polynomial of degree 
two in x 2 : Q{x) = x 2M+2 R(x) where R(x) = (qM+3 x A + qjyi+2 x 2 + qu+i)- 
This can be done if we choose the p's in order to cancel all the q m for m ^ M. 
For m G {0, . . . , M — 1}, equations ([TBI) determine uniquely all the p's up to 
a common factor po than can be taken to one without loss of generality: p n 
is a polynomial in (E, A) of degree at most n. From equation (|16|) written 
for m = M, the condition qu = imposes an algebraic relation of degree 
at most (M + 1) between E and A that implicitly defines A as a function 
of E (several branches are possible in general). After eliminating A, the 
discriminant of R appears as a function Am(E,s) of E and s only. For a 
given s, all the values of E where Am < are forbidden. The boundaries 
of the gaps are given by the zeroes of E t— ► Am(E,s). For instance, when 
M = we have: 

A (E, s) = 65536 E 6 - 73728 E 4 s - 41472 E 3 + 20736 E 2 s 2 + 7776 Es + 6561. 

(17) 

For a given value of s, the real roots of this polynomial define intervals of 
forbidden values of E. The degree of the polynomial Am(E,s) increases 
with M. 



n 





1 


2 


3 


4 


5 


6 


E n 


-0.5012 


-0.2549 


0.9606 


2.1003 


3.5281 


5.1202 


6.8609 



M 


E < ■ ■ ■ < Ei 


Ei < • • • < Ei 


E 2 < ■■■ < E 3 


E 3 < ■■■ < E 4 


E4. < ■ ■ ■ < E 5 


E 5 < ■ ■ ■ <E 6 


<4 
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[0.0399,0.4582] 










6 




[-0.1108,0.4359] 
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[-0.1836,0.3908] 


[1.1734,1.4748] 
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[-0.4614,-0.4281] 


[-0.2001,0.3408] 


[1.0601,1.5156] 








15 


[-0.4994, -0.4865] 
[-0.4835, -0.4336] 


[-0.2414,0.0507] 


[0.9845,1.2930] 


[2.1214,2.3143] 
[2.3789,2.7014] 


[3.7509,4.3535] 


[5.3268,6.0726] 


16 


[-0.5006, -0.4369] 


[-0.247,-0.24] 
[-0.2325,0.0206] 


[0.9903,1.2652] 


[2.1318,2.3167] 
[2.3307,2.6590] 


[3.6068,3.7792] 
[3.9847,4.3128] 


[5.3426,6.0835] 


16 

u 

M=0 


[-0.5006,-0.4211] 


[-0.2497,0.4582] 


[0.9726,1.5156] 


[2.1214,2.8710] 


[3.6068, 4.4028] 


[5.3268,6.0835] 



Table 1: Comparison between the exact energies {E n \ n = 0, 1, 2, . . . } (given 
in the upper table) and the gaps for s = —2.3 and several values of M. 
The two lowest energies (Eq,E\), being negative though below the central 
barrier, form a tunnelling doublet. There are two gaps between (Eq,Ei) 
for M = 15, and one gap for M = 16. The last line sums up the best 
bounds when comparing the gaps obtained for each M up to 16. We did 
not retain the lowest bounds that are below the trivial bound miny. 

We can generalize and reformulate these algebraic manipulations to any 
kind of polynomial potential of even degree degV (not necessarily sym- 



6 



-6 -5 -4 -3 -2 -1 1 2 3 4 5 6 -6 -5 -4 -3 -2 -1 1 2 3 4 5 6 
S S 




Figure 1: (color on line) The black thick lines stand for the exact spectrum 
of the Hamiltonian p 2 /2 + sx 2 /2 + x 4 /4. On the upper left pannel the 
boundaries of the gaps that are given by the solutions of (I18p for M = (red, 
the equation of the curve is obtained by cancelling flTTJ)), M = 1 (green), 
M = 2 (blue). The gaps on the upper right pannel are obtained for M = 7 
(magenta) and M = 15 (cyan). For a given M, no gap can be found for a 
too small value of s but we can see that for a given value of s all the lowest 
energies can be separated by a gap for a large enough M. The boundaries 
of the gaps converge to the exact spectrum when M increases. For instance, 
for M = 15 and s > 0, the lowest gap is very thin and can hardly be 
distinguished from the ground state Eq. In some cases, the equations (fT8|) 
provide two disconnected gaps in between two successive energies. As can be 
seen in the pannel below, for M = 7 and M = 15, even the doublet (Eq, E\) 
can be separated in a tunnelling regime (i.e. when they are both negative). 

metric). With the choice Q12|). the expression (jSJ) shows that i^e^ 2 / 2 
is a polynomial Q(x,X,E) of degree degV + degP + 1 in x. For any 
(A, E) the coefficients of P can be chosen to cancel all the coefficients 
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of x i ^ Q(x,X,E) but the (degV + 2)th higher powers. The cancella- 
tion of the coefficient q(X,E) of x degP allows the factorization Q(x, \,E) = 
x degP+1 R(x, A, E) where the degree of x i— ► R(x,X,E) is precisely degF. 
Now, if degP is odd and x i— ► R(x,X,E) has no root of odd multiplicity 
for a given couple (A, E) such that q(X, E) = 0, we are sure that E can- 
not be an eigenvalue. In the (x, X, £')-space, the two equations q(X, E) = 
0, R(x, X, E) = define generically a curve C (possibly made of disjoint 
smooth pieces) that can be parameterized by x, namely (E{x), A(x)), where 
the Jacobian J = \d\q&ER — d\RdEq\ does not vanish. The projection 
of C on the E-axis defines some intervals outside which no eigenenergy can 
be found. The borders of these intervals are the projections of some points 
(not necessarily of all points) of C where the tangent is normal to the E-axis. 
Using the implicit function theorem these points are to be found among the 
solutions of = dE/dx = —d\qd x R/ J. To sum up, the boundaries of the 
gaps are to be found from the solutions — if there are any — of the three 
equations 

q(X,E) = 0, R(x,X,E)=0, d\q(X, E) d x R(x, A, E) = 0. (18) 

The boundaries obtained above for the quartic potential using the zeroes 
of the discriminants Am are included in the solutions of (|18|) since they 
cancel simultaneously q, R and d x R but some others boundaries may be 
obtained if d\ q vanishes instead of d x R. Determining which values represent 
boundaries of gaps from all the solutions of (|18p may require some global 
analysis that can be pursued numerically. The larger M, the larger the 
degree of the algebraic equations to solve and the more gaps are expected 
to be found. For the quartic oscillator, the numerical results are shown for 
several values of M in figure [1] and table Q] illustrates the precision of our 
method. 

Let us consider also the even potential V(x) = x 6 . The polynomial R 
is now of third degree in the variable x 2 . As explained above, a necessary 
condition to obtain the boundaries of the gaps for a given positive integer 
M is to solve simultaneously the three equations (fTBj) . The last step consists 
in verifying that the polynomial R(x) does not have real roots of odd multi- 
plicity. The results obtained are qualitatively similar that the one's showed 
above for the quartic oscillator (see tables [2]) . 

4 Conclusions 

From what precedes it is clear that our method is not specific to the quartic 
potential but can be adapted straightforwardly to many other situations. 

This method can be also applied to a spectral problem with different 
boundary conditions. For, say, Dirichlet boundary conditions equation (ii) 
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n 





1 


2 


E„ 


0.6807 


2.5797 


5.3948 



M 


Eq < ■ ■ ■ < i?i 


E\ < ■ - ■ < E2 









1 




0.967,1.041 






2 




0.736; 1.122 






3 




0.710,1.081 




[3.016,3.624] 


4 




0.715,1.036 




[2.724,3.702 


5 




0.723,0.995 




[2.674,3.642] 



Table 2: Comparison between the exact energies {E n \ n 
the left table) and the gaps for V(x) 



x 6 . 



0, 1, 2} (given in 



is replaced by the vanishing of Jn at the two points where the function (p 
vanishes. 

Even though our method cannot provide a rigorous proof of the exis- 
tence of an eigenenergy in a given range, it offers some precise clues where 
some possible energies may lie, specially if some convergence behaviour is 
observed, as it appears in the case of the quartic oscillator. Therefore we 
are able to obtain this way a substantial piece of physical information. 

Aknowlegments: We would like to thank Jeremy Ledeunff for helping 
us to check the eigenenergies for the x 6 -potential. 
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